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Abstract 

A  comparison  of  the  performance  of  a  solid-polymer  electrolyte  direct  methanol  fuel  cell  (SPE-DMFC)  with  aqueous  methanol  and  air 
mixed  at  the  anode,  and  with  only  aqueous  methanol  is  presented.  The  performance  of  the  mixed  methanol  reactant  and  air  fed  anode  SPE- 
DMFC  is  superior  to  the  conventional  methanol  only  fed  anode  SPE-DMFC.  The  experimental  performance  of  the  mixed-reactant  anode  SPE- 
DMFC  are  compared  with  results  of  a  numerical  modelling  based  on  phenomenological  transport  equations  for  electrocatalyst  layer,  diffusion 
layer  and  the  polymer  electrolyte  membrane. 

©  2002  Published  by  Elsevier  Science  B.V. 
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1.  Introduction 

The  direct  use  of  methanol  in  fuel  cells  is  considerably 
attractive  from  the  point  of  view,  of  simplicity  of  system 
design  and  hence  cost  [1].  A  direct  methanol  fuel  cell 
(DMFC)  consists  of  an  anode  at  which  methanol  is  electro- 
chemically  oxidised  to  carbon  dioxide  according  to 

CH3OH  +  H20  ->  C02  t  +  6H+  +  6e“  (1) 

and  a  cathode  at  which  oxygen  is  reduced  to  water 

1 02  +  6H+  +  6e“  ->  3H20  (2) 

Accordingly,  the  overall  cell  reaction  is 

CH3OH  +  §02^C02T  +  2H20  (3) 

The  present  day  DMFCs  employ  a  thin  proton-conducting 
polymer  sheet  as  the  electrolyte  and  hence  are  commonly 
referred  to  as  solid  polymer  electrolyte  direct  methanol  fuel 
cells  (SPE-DMFCs).  During  the  last  decade,  significant 
advances  have  been  made  in  the  DMFC  development.  Power 
densities  of  450  and  300  mW/cm2  under  oxygen  and  air  feed 
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operation,  respectively,  and  200  mW/cm2  at  a  cell  potential 
of  0.5  V  have  been  reported  for  cells  operating  at  tempera¬ 
tures  close  to  or  above  100  °C  under  pressurised  conditions 
with  platinum  loadings  of  1-2  mg/cm2  [2,3].  Besides,  the 
development  of  DMFC  stacks  for  both  transportation  and 
portable,  applications  has  gained  momentum  in  the  last  2-3 
years,  and  stacks  with  power  densities  of  1  kW/dm3  and 
an  overall  efficiency  of  37%  at  a  design  point  of  0.5  V 
per  cell  have  been  accomplished  [4].  The  performance  of 
SPE-DMFCs  is  thus  competitive  with  respect  to  the  reformer- 
based  H2/air  polymer  electrolyte  fuel  cells  (PEFCs),  espe¬ 
cially  if  one  considers  the  complexity  of  the  latter  whole 
system  [5].  In  a  recent  review  on  the  prospects  of  fuel  cells 
for  road-transportation  purposes  [6],  it  has  been  suggested 
that  only  fuel  cell  vehicles  (FCVs)  with  directly  fuelled  fuel 
cells,  such  as  the  SPE-DMFCs,  will  be  able  to  exceed  the 
performance  of  future  conventional  internal  combustion 
engine  vechiles  (ICEVs).  However,  further  improvements 
in  the  anode  performance  of  SPE-DMFCs  in  addition  to  a 
reduction  in  their  cost  would  be  mandatory  for  their  use  in 
FCVs  [7].  A  step  in  this  direction  appears  to  be  the  devel¬ 
opment  of  mixed-reactant  SPE-DMFCs  [8,9].  It  is  note¬ 
worthy  that  selective-electrode  fuel  cells  with  mixed- 
reactant  feed  have  been  considered  for  solid  oxide  [10] 
and  proton-exchange  membrane  fuel  cells  [11,12]. 

In  this  communication,  we  report  the  performance  of  a 
SPE-DMFC  with  a  mixed-reactant  anode  in  conjunction 
with  its  numerical  modelling.  The  study  demonstrates  that 
the  SPE-DMFC  with  mixed-reactant  anode  exhibits  superior 
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Nomenclature 

a  effective  catalyst  area  per  unit  volume 

Ad  anode  diffusion  region 

AR  anode  reactive  region 

ct  concentration  of  species  i 

Co  cathode  diffusion  region 

CR  cathode  reactive  region 

Dt  diffusion  coefficient  of  species  i 
Dt_j  pair  diffusion  coefficient  of  gas  species  i  and  j 

f  Faraday  constant  in  units  of  RT 

F  Faraday  constant 

/i  ionic  current  density 

4  electronic  current  density 

*o,ref  exchange  current  density  at  reference  condi¬ 

tions 

I  operative  cell  current  density 

j  transfer  current  per  unit  volume 

kp  hydraulic  permeability 

electrokinetic  permeability 
Ki  Henry’s  constant  for  species  i 

M  membrane  region 

n  number  of  electrons 

Ni  mass  flow  of  species  i 

P  pressure 

R  gas  constant 

Si  stoichiometric  coefficient  of  species  i 

S\  liquid  saturation 

T  temperature 

v  water  velocity 

Xi  molar  fraction  of  species  i 

z  fixed-site  charge 

Greek  symbols 

aa  anodic  transfer  coefficient 

ac  cathodic  transfer  coefficient 

em  membrane  porosity 

£p  diffusion  region  porosity 

<pi  electric  potential  of  liquid  phase 
(j)s  electric  potential  of  solid  phase 

y  kinetic  factor 

v\  interface  overpotential 

k  ionic  conductivity 

/j  water  viscosity 

p  water  density 

a  electronic  conductivity 

Subscripts  and  superscripts 
a  anode 

av  average 

c  cathode 

eff  effective  value 

f  fixed-charge 

g  gas  phase 

in  inlet 


1  liquid  phase 

out  outlet 

ref  reference  conditions 

sat  saturated  vapor 

w  water 


performance  in  relation  to  conventional  selective-anode 
SPE-DMFCs. 


2.  Experimental 

A  SPE-DMFC  was  assembled  with  a  Nafion®-117  based 
membrane  electrode  assembly  (ME A)  obtained  from  Johnson 
Matthey,  which  employed  Pt-Ru/C  as  the  anode  catalyst  with  a 
loading  of  1  mg/cm2  of  Pt  and  0.46  mg/cm2  of  Ru,  and  Pt/C 
as  the  cathode  catalyst  with  a  loading  of  4.6  mg/cm2  of  Pt.  In 
the  experimental  cell,  the  anode  and  cathode  of  the  MEA 
were  in  contact  with  high-density  graphite  blocks,  impreg¬ 
nated  with  phenolic  resin,  into  which  were  cut  the  gas/liquid 
flow  channels  [13].  The  ridges  between  the  channels  pro¬ 
vided  the  electrical  contact  to  the  carbon  electrodes  and  the 
total  machined  geometrical  area  of  25  cm2  was  taken  as  the 
active  area  of  the  cell.  Electrical  heaters  were  mounted  at  the 
rear  of  the  graphite  blocks  to  maintain  the  desired  cell 
temperature  which  was  controlled  through  a  temperature 
controller  and  monitored  by  thermocouples  buried  in  the 
graphite  blocks.  Gold-plated  metallic  bolts  were  screwed 
into  the  blocks  to  allow  electrical  contact.  For  selective- 
anode  tests,  a  solution  of  1 .0  mol/dm3  (M)  aqueous  metha¬ 
nol  was  fed  to  the  anode  at  a  feed  rate  of  ca.  15  cm3/min 
and  the  unreacted  fuel  returned  to  the  storage  reservoir 
where  the  C02  gas  product  was  released  to  the  atmosphere. 
For  the  mixed-reactant  anode  test,  either  nitrogen  gas  or  air 
pre -heated  to  cell  temperature  was  fed  at  200  cm3/min  to  the 
anode  along  with  1  M  aqueous  methanol. 

On  polarising  the  anode  galvanostatically,  hydrogen  was 
evolved  in  accordance  to  the  reaction: 

6H+  +  6e_  ->  3H2  f  (4) 

at  the  cathode  which  also  served  as  a  reference  electrode 
[8,14].  To  obtain  the  cell  polarisation  data,  hot  air  at 
200  cm3/min  was  also  supplied  to  the  cathode  from  an  air 
cylinder.  The  experimental  set-up  employed  for  the  study  is 
shown  schematically  in  Fig.  1. 

3.  Model  description 

The  model  uses  a  cell  configuration  comprising  five 
regions  as  shown  in  Fig.  2,  each  of  which  has  a  different 
chemical  and  physical  transport  phenomenon.  The  model  is 
primarily  a  modification  of  the  approach  described  earlier 
for  modelling  a  polymer  electrolyte  fuel  cell  [15-18].  In 
brief,  the  model  uses  a  one-dimensional  approach  where  the 
independent  variables  are  the  electrical  potential  of  the  solid 
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Fig.  1.  Schematic  diagram  of  the  experimental  set-up. 


(f/)s)  and  liquid  (</> i)  phases,  the  hydraulic  pressure  (P), 
concentrations  of  methanol  (c'ch3oh),  carbon  dioxide 
(cco2),  and  oxygen  (c02)  in  the  liquid  phase,  and  nitrogen 
molar  fraction  (jcn2)  in  the  gas  phase.  The  modelled  transport 
phenomena  and  conservation  equations  are  in  given  in 
Table  1.  Any  spatial  variation  in  the  variables  is  governed 
by  the  respective  phenomenological  expressions,  which 
relate  the  flux  of  any  physical  quantity  to  its  generating 
force.  The  closed  set  of  differential  equations  for  the  model, 
obtained  by  combining  the  phenomenological  and  conser¬ 
vation  equations,  are  described  in  Table  2. 


The  model  uses  (a)  the  Stefan-Maxwell  equation  to 
account  for  gas-diffusion  in  porous  medium,  (b)  the 
Nemst-Planck  equation  to  account  for  the  proton  conduc¬ 
tion  and  transport  of  the  neutral  chemical  species  in  the 
electrolyte  and  the  bulk-liquid  phase,  (c)  Ohm’s  law  for 
electronic  conductivity  in  the  bulk- solid  phase,  (d)  a  mod¬ 
ified  Schlogl  equation  [19,20]  to  account  for  convective 
motion  in  the  membrane  electrolyte  and  ME  A  comprising 
the  backing  plus  diffusion  layers  where  it  takes  the  form  of 
the  Darcy  equation.  The  model  source  equation  for  the 
reaction  kinetics  at  the  anode  and  cathode  are  obtained 


Regions 


I  II  III  iv  v 
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Fig.  2.  Cell  configuration  adopted  for  numerical  modelling. 
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Table  1 

Model  transport  phenomena  and  its  conservation  equations 

Regions 

Variables 

Flux 

Conservation 

CCH3OH,  CC02,4>s:  ? 

VcH3OH:  VC02, 4,  V 

V  •  AZch3oh  =  V  •  Nqo2  =  V  •  4  =  V  •  v  =  0 

CCH3OH,CCO2,0S,01,E 

Vch30H;  Vco2, 4,  i\,v 

V  •  Vch,oh  =  ycH30H,  V  •  NCo2  =  yco2 ,  V  •  4  =  - V  •  i\  =  /,  V  •  v  =  pymo 

M 

CCH3OH,Co2,  01,  P 

Vch3oh- Vo2,  zi,  y 

V  •  AZch3oh  =  V  •  No2  =  V  •  i\  =  V  •  v  =  0 

CR 

^ch3oh,  Od2  j  0S,  0i,  p 

Vch3oh- V02, 4,  z'i,  v 

V  ■  /Vo,  =  y02 ,  V  ■  4  =  -V  •  h  =/,  V  •  v  =  pv’h2o 

Cd 

Xn2,Xo2,0s  1P 

Vn2,Vo2, 4,  v 

V  •  Vn2  =  V  •  Nq2  =  V  •  4  =  0,  V  •  v  oc  V  •  ZVh2o 

In  the  table  y,  =  (sf  /n^F)f,  where  (3  =  a,c. 


electrolyte  membrane,  (e)  variable  gas-liquid  saturation  on 
the  cathode,  and  (f)  liquid-phase  transport  at  the  anode.  Since 
the  catalyst  loading  at  the  cathode  is  substantially  high, 
any  methanol  crossover  effect  on  the  cathode  polarisation 
has  been  neglected  in  the  numerical  modelling  of  the  cell. 

Since  the  present  model  is  based  on  a  mechanistic 
approach,  it  requires  knowledge  of  several  material  pro¬ 
perties,  e.g.  effective  diffusion  coefficients,  hydraulic  and 
electro-kinetic  permeability,  proton  conductivity  in  the  elec¬ 
trolyte,  electronic  conductivity  in  the  electrodes,  etc.  These 
experimental  parameters  have  been  taken  from  the  literature. 
Only  electronic  conductivity  and  the  pre-exponential  kinetic 
factor  have  been  appropriately  adjusted  in  the  numerical 
modelling.  Besides,  the  estimated  effective  values  of  all 
the  transport  coefficients  account  for  the  porosity  of  the 
medium  under  consideration.  The  generic  effective  diffusion 
and  permeability  coefficients,  keff  and  Z)eff,  respectively,  for 
various  reactant/product  species  in  the  membrane  electrolyte 
are  obtained  by, 

keff  =  k8m  (7) 


Variables 

Equations 

Validity  Regions 

Phenomenological  equation 

cch3oh 

V  •  (£>ch3ohVcch3oh  -  cch3ohv)  =  0 

M 

Nemst-Planck 

V  •  (£>ch3ohVcch3oh  —  cch3ohv)  =  Very  oh 

An 

cco2 

V  •  (£>co2VcCo2  -  cco2v)  =  0 

M 

Nemst-Planck 

V  •  (Dco,  VcCo2  -  cco2v)  =  >’co2 

An 

co2 

V  •  (Do,  Vco2  —  co2v)  =  0 

M 

Nemst-Planck 

V  •  (D0,Vco2  -  c02v)  =  yCo2 

Cn 

Xn2 

VxN2  =  ^;i(RT/PDN2_,)(xN2AZ,-x,AZN2) 

where  j  =  02,  H20 

Cn 

Stefan-Maxwell 

*02 

VxG2  =  J2nj=l(RT/PD02-j)(x°2Nj-xjNo2) 

where  j  =  N2,  H20 

Cn 

Stefan-Maxwell 

0s 

V  •  (<rV0s)  =  0 

M 

Ohm’s  law 

V  •  (uV0s)  =/ 

CA 

V  •  (<rV0s)  =/ 

Cn 

01 

V  •  (kV0!  +  FCfEv)  =  0 

M 

Nemst-Planck 

V  •  (kV0!  +  FcfEv)  =f  -  (Fcf/p)yu2o 

CA 

V  •  (kV0!  +  FcfSv)  =f  -  (Fcf/p)yR2o 

Cn 

p 

V  •  (~(kj,/p)VP)  =  0 

CD 

Darcy 

V  •  (6V0!  -  (kF/p)VP)  =  0 

M 

Modified  Schlogl 

V  •  (ZA70!  -  (kP/p)VP)  =  -(l/p)yu2o 

CA 

V  •  (W0!  -  (kP/p)S7P)  =  -(1/pW 

Cn 

using  Butler-Volmer  equations,  which  assume  the  Tafel 
form  for  irreversible  reactions,  as  described  below. 


V  •  z a  =  aio  s 


cch3oh 


(5) 


where 

/  =  -*; 

J  RT, 

and 


*7a  =  0s, a  -  01, a 


(6) 


V.ic  =  aio,c^)  [-e^H 

where 

=  0S, c  -  01, c 


The  other  characteristics  of  the  model  assume  (a) 
steady-state  operation  of  the  cell,  (b)  constant  temperature 
(isothermal)  operation  of  the  cell,  (c)  zero  differential 
gas  pressure  in  the  porous  medium,  (d)  fully-hydrated 

Table  2 

System  of  coupled  differential  equations  for  the  numerical  modelling 


In  the  table  b  =  {k^J p)zfC{F . 
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and 

Deff  =  Dem  (8) 

where  em  is  the  membrane  fraction  in  the  catalyst  layer. 

In  the  porous  region,  the  total  porosity  (ep)  may  be 
partially  occupied  by  water  and  gas  and  therefore,  we  need 
to  specify  the  relative  occupancy  of  either  phase  while 
accounting  for  the  effective  diffusion  and  permeability 
coefficients  at  the  anode  and  cathode.  In  the  electrode  porous 
region,  the  Bruggemann  correction  to  the  diffusion  coeffi¬ 
cients  is  applied  for  porous  systems  with  homogeneous 
porosity  and  particle  size  [21]. 


Accordingly,  at  the  anode, 

D'“  =  D(4)13  (9) 

where  ej,  =  spSi,  Si  being  the  liquid-phase  saturation  para¬ 
meter.  Similarly,  at  the  cathode, 

D‘tt  =  D(£«)15  (10) 

where  s®  =  £p(l  —  Si).  The  hydraulic  permeability  (kp) 
accounts  for  changes  in  liquid  saturation  in  accordance  with 
the  equation  [22,23] 

*f  =  M?  (ii) 

The  full  set  of  model  parameters  is  given  in  Table  3. 


Table  3 

Model  parameters  employed  in  the  numerical  modelling 


Parameters 

Values 

Reference 

Physical  and  operative  parameters 

Membrane  thickness 

230  pm 

Diffusive  region  thickness 

260  pm 

Catalyst  layer  thickness 

10  pm 

Cell  temperature  (7) 

80  °C 

Anode  inlet  pressure  P'f 

1  atm 

Cathode  inlet  pressure  P[n 

1  atm 

Nitrogen-oxygen  mole  ratio 

0.79/0.21 

Anode  feed  flow  rate  (liquid) 

15  ml/min 

Cathode  feed  flow  rate  (gas) 

200  ml/min 

Membrane  parameters 

Ionic  conductivity  ( k ) 

0.14  mho/cm 

[25] 

CH3OH  diffusivity  (£>Ch3oh) 

2.59  x  10~6  cm2/s 

[25] 

C02  diffusivity  (DCXh  ) 

2.5  x  10-6  cm2/s 

02  diffusivity  (D0,) 

1.22  x  10~6  cm2/s 

[26] 

Fixed  charge  concentration  (cf) 

1.2  x  10  3  moFcm3 

[18] 

O2  Henry’s  constant  (Kq2  ) 

2.02  x  105  atm  cm3/mol 

[26] 

Hydraulic  permeability  (kp) 

1.8  x  10~14  cm2 

[18] 

Electrokinetic  permeability  (kf) 

7.18  x  10~16  cm2 

[18] 

Water  viscosity  ( p ) 

4.56  x  10-4  kg/m  s 

[27] 

Water  density  ( p ) 

0.054  mol/cm3 

[27] 

Electrode  parameters 

Electronic  conductivity  ( 0 ) 

0.38  mho/cm 

CH3OH  diffusivity  (in  water,  Dch,oh) 

8.5  x  10  5  cm2/s 

[28] 

C02  diffusivity  (in  water,  Dqq2) 

3.5  x  lO-5  cm2/s 

O2-N2  pressure-diffusivity  (Do2-n2) 

0.279  atm  cm2/s 

[18] 

02-H20  pressure-diffusivity  (Dq2_h2o) 

0.370  atm  cm2/s 

[18] 

H20-N2  pressure-diffusivity  (Dh2o-n2) 

0.387  atm  cm2/s 

[18] 

Membrane  fraction  in  catalyst  layer  (em) 

0.4 

[18] 

Diffuser  layer  porosity  (ep) 

0.4 

[18] 

Liquid  saturation  in  anode  (S£) 

0.75  mixed-reactants  feed,  0.50  selective-reactant  feed 

Electrokinetic  permeability  (kf) 

7.18  x  10-16  cm2 

[18] 

Kinetics  parameters 

Pre-exponential  parameter  (m'o.rcf) 

2.6  x  10  A/cm3  anode,  2.0  x  1 0  4  A/cm3  cathode 

Cathodic  transfer  coefficient  (ac) 

1.2  for  oxygen  reduction 

[16] 

Anodic  transfer  coefficient  (aa) 

0.5  for  methanol  oxidation 

[29] 

02  reference  concentration  (co2,ref) 

4.72  x  10  7  mol/cm3 

CH3OH  concentration  (ccH3OH,ref) 

1.0  x  10-3  moFcm3 

02  concentration  parameter  (y02 ) 

1.0 

[16] 

CH3OH  concentration  parameter  (/ch3oh) 

1.0 

[29] 

CH3OH  stoichiometric  coefficient  (.s’ch3oh) 

1 

C02  stoichiometric  coefficient  (sC(h ) 

-1 

H20  stoichiometric  coefficient  G?h2o) 

1  anode,  3  cathode 

O2  stoichiometric  coefficient  (,s’o2 ) 

-3/2 

Number  of  electrons  ( n ) 

6 
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Table  4 

Internal  boundary  conditions  for  the  system  of  coupled  differential  equations 


AdMr 

ar/m 

M/Cr 

cR/cD 

CC.H3OH 

A;h3ohVcch3oh  -  = 

^ch3oh  V  CcH3OH  + 

A:h3ohVcch3oh  _  = 

A:h3  oh  VcCh3oh|+ 

A:h3ohVcch3oh  _  = 
•Dch3ohVcch3oh|+ 

cch3oh  -  =  0 

cco2 

Dc02Vcco2  -  = 

A;o2  Vcco  + 

cco2  -  =  0 

~ 

co2 

XN2a 
*02  3 

co2  +  =  0 

£>o2Vc02|_  =  Do,  Vco2  + 

co2|+  =xo2\+Pm/Ko2 

- 

- 

- 

- 

</>s 

<tV</>s  =  + 

V0S|_  =0 

V<K|+  =o 

<A7(/>S  _  =  crV(/>s  + 

0il+ 

kV</>!  _  =  kV^1  + 

kV^|_  =  KV(f>i\+ 
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In  the  table  b  =  / /j.)z(Ci  F:  the  subscripts  —  and  +  refer  to  the  regions  on  the  left  and  on  the  right  of  the  interface,  respectively. 

a  The  boundary  conditions  are  given  at  the  CD/CF  interface. 


Solution  of  the  numerical  model  is  obtained  once  a  proper 
set  of  boundary  condition  is  given.  For  one-dimensional 
models,  where  cell  inlet  and  outlet  are  coincident,  the  task  is 
not  trivial  and  a  number  of  assumptions,  valid  for  the  system 
under  consideration,  are  desired.  For  example,  in  the  cell,  we 
assume  that  hydraulic  inlet  pressure  equals  its  outlet  pres¬ 
sure,  i.e.  Pm  =  Pout.  Besides,  for  the  chemical  species  con¬ 
centration,  it  is  mandatory  to  obtain  a  mass  balance  between 
the  inlet  and  outlet  compositions  of  the  reactants/products. 
While  the  inlet  compositions  are  known,  the  outlet  composi¬ 
tions  are  unknown  a  priory.  Therefore, 

cav  =  i  (cf  +  c°ut)  (12) 

has  to  be  estimated  by  an  iterative  procedure  which  yields 
self-consistent  values  for  c°ut  and  its  flux  A°ut.  For  the 
liquid-phase  transport  at  the  anode  side, 
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(13) 


where  cw  a  is  the  water  concentration  assumed  to  be  constant 
as  water  is  the  excess  solvent,  Nff  and  A°u*  are  the  outlet 
flux  of  the  species  i  and  water,  respectively.  Similarly,  for  the 
gas-phase  transport  at  the  cathode  side,  the  molar  fraction  of 
the  species  at  the  outlet  is, 
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(14) 


Further  details  on  mass  balance  for  fuel  cell  systems  can 
be  found  elsewhere  [24].  The  cathode  solid-phase  potential 
in  the  catalyst-diffuser  region  interface  is  obtained  by 
using 


0s  =  01 


1 

&cf 
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derived  by  an  approximate  analytical  integration  of  the 
Butler- Volmer  equation  over  the  entire  catalyst  layer  volume 
[15].  The  approximations  applied  for  the  Butler- Volmer 
equation  integration  at  the  cathode  side  are  not  valid  at  the 
anode.  Therefore,  we  are  not  able  to  write  an  explicit  relation 
for  the  potential  at  the  catalyst-diffuser  region  interface  and 
an  iterative  optimisation  procedure  is  required  to  obtain  its 
value.  For  the  liquid-phase  potential,  an  arbitrary  value  can  be 
given  at  the  anode  catalyst-diffuser  region  interface.  The 
internal  boundary  conditions  employed  to  couple  the  five 
numerical  regions  and  to  solve  the  entire  set  of  equations  are 
given  in  Table  4. 


4.  Results  and  discussion 

The  galvanostatic  polarisation  data  for  the  selective  and 
mixed-reactants  anode  tests  are  given  in  Fig.  3.  Unlike  the 
data  reported  by  Barton  et  al.  [8]  for  their  mixed-reactant, 
strip-cell  DMFC,  we  find  a  substantial  improvement  in 
the  performance  of  the  anode  particularly  at  load  current 
densities  beyond  ~400  mA/cm2.  But,  akin  to  the  findings 
of  Barton  et  al.  [8],  we  also  observe  little  difference 
between  the  anode  polarisation  data  obtained  with  the 
mixed  feed  of  1  M  aqueous  methanol  plus  air,  and  the 
anode  polarisation  data  obtained  with  1  M  aqueous  metha¬ 
nol  plus  nitrogen.  These  data  clearly  reflect  that  the  added 
oxygen  does  not  lead  to  any  significant  parasitic-oxidation 
of  methanol  at  the  anode.  Besides,  unlike  the  observation 
of  Barton  et  al.  [8],  we  do  not  find  the  air  stream  to  impede 
the  mass  transfer  of  methanol  to  the  anode  surface.  Instead, 
we  observe  that  air-feed  helps  scavenge  carbon  dioxide 
from  the  active  catalytic  sites  ameliorating  oxidation  of 
methanol  on  the  anode. 

It  is  worth  noting  that  a  SPE-DMFC  operating  at  1  A  load 
would  require  5.8  x  10-5  dm3/s  of  oxygen  at  the  cathode, 
and  7.06  x  10-8  dm3/s  of  methanol  at  the  anode,  and  will 
result  in  water  production  of  6.26  x  10-8  dm3/s  at  the 
cathode  and  C02  exhaust  of  3.87  x  10-5  dm3/s  at  its  anode. 
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Fig.  3.  Anode  polarisation  data  obtained  at  80  °C  for  (■)  selective  anode  with  1  M  aqueous  methanol  at  feed  rate  of  ca.  15  cm3/min,  (#)  mixed-reactant 
anode  with  1  M  methanol  at  feed  rate  of  ca.  15  cm3/min  plus  200  ml/min  air  at  atmospheric  pressure,  and  (A)  mixed-reactant  anode  with  1  M  methanol  at 
feed  rate  of  ca.  15  cm3/min  plus  200  cm3/min  nitrogen  at  atmospheric  pressure. 


This  represents  ca.  550-fold  volume  increase  in  the  anode 
compartment  of  the  cell  during  its  operation  and  suggests 
that  efficient  C02  removal  from  the  anode  would  be  seminal 
to  further  the  cell  performance. 

Galvanostatic  polarisation  data  for  the  complete  cell  are 
given  in  Fig.  4.  A  high  catalyst  loading  of  4.8  mg/cm2  of  Pt 
on  the  cathode  was  used  intentionally  to  minimise  the 
cathode  losses  owing  to  any  methanol  crossover  from  the 
anode  to  the  cathode.  The  data  clearly  suggest  that  the  cell 
performance  with  mixed-feed  anode  to  be  superior  in  rela¬ 
tion  to  its  polarisation  data  for  the  methanol  only  fed  anode. 


The  polarisation  plots  obtained  from  the  numerical  mod¬ 
elling  of  the  SPE-DMFC,  its  anode,  both  with  selective  and 
mixed-reactant  configurations,  and  its  cathode  are  given  in 
Figs.  5-7.  A  comparison  of  the  anode  polarisation  data 
obtained  from  the  experiments  (Fig.  3)  and  modelling 
(Fig.  6)  suggests  them  to  be  in  conformity.  The  only  para¬ 
meter  that  differs  in  the  numerical  modelling  of  the  polar¬ 
isation  data  for  the  selective  and  mixed-reactant  anodes  is 
the  liquid  saturation  (Si)  in  the  diffusion  layer  of  the  anode 
which  has  the  relative  values  of  0.5  and  0.75,  respectively.  In 
the  mixed-reactant  anode,  carbon  dioxide  evolving  at  the 


Fig.  4.  Cell  polarisation  data  obtained  at  80  °C  for  (■)  selective  anode  with  1  M  aqueous  methanol  at  feed  rate  of  ca.  15  cm3/min  and  cathode  air  flow  rate  of 
200  cm3/min  at  atmospheric  pressure,  and  (A)  mixed-reactant  anode  with  1  M  methanol  at  feed  rate  of  ca.  15  cm3/min  plus  200  cm3/min  air  at  atmospheric 
pressure  and  cathode  air  feed  at  200  cm3/min  at  atmospheric  pressure. 


50 


A.K.  Shukla  et  al./ Journal  of  Power  Sources  111  (2002)  43-51 


Current  density  /  mAcm'2 

Fig.  5.  Numerical  modelling  for  polarisation  of  the  SPE-DMFC  with  selective  and  mixed-reactant  anodes  akin  to  the  experimental  conditions. 


anode  will  be  easily  removed  from  the  electrode  pores,  as  it 
will  readily  mix  with  the  air  stream.  By  contrast,  in  the 
selective -reactant  anode,  evolving  carbon  dioxide  will  reg¬ 
ister  a  back-pressure  in  the  electrode  pores  due  to  the  liquid 
stream  in  the  flow  field  which  will  hinder  its  escape  from  the 
pores.  Thus,  a  reduction  in  the  liquid- saturation  parameter 
will  obstruct  the  effective  supply  of  methanol  to  the  elec¬ 
trode  surface  bringing  about  its  polarisation  much  earlier  in 
relation  to  the  mixed-reactant  anode. 


From  the  numerical  modelling  of  the  air  cathode  (Fig.  7), 
the  presence  of  an  excessive  amount  of  saturated  water 
vapour  (x^at)  in  electrode  pores  in  the  air  cathode  at  80  °C 
and  1  atm  of  air  seems  to  be  the  limiting  factor  for 
its  performance.  This  suggests  that  an  increase  in  either 
the  flow  rate  of  the  air  feed  at  the  cathode  or  its  pre¬ 
ssure  will  be  mandatory  to  further  the  performance  of  the 
SPE-DMFC.  There  is  no  reason  to  believe  that  the  cell 
engineering  and  its  operating  conditions  are  yet  fully 


Fig.  6.  Numerical  modelling  of  the  anode  polarisation  for  selective  and  mixed-reactant  anodes  akin  to  the  experimental  conditions. 
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Fig.  7.  Numerical  modelling  of  the  air-cathode  polarisation  at  air-feed  rate  of  200  cm3/min  at  atmospheric  pressure  akin  to  the  experimental  conditions. 


optimised,  and  further  improvements  in  the  cell  perfor¬ 
mance  are  highly  likely. 


5.  Conclusions 

This  study  has  shown  that  using  a  mixed  air  and  methanol 
solution  anode  stream  in  the  DMFC  leads  to  an  increase  in 
anode  and  cell  performance.  It  is  suggested  that  this 
improvement  is  due  to  a  higher  liquid  saturation  in  the 
anode  diffusion  layer  of  the  DMFC  and  the  faster  removal 
of  carbon  dioxide  at  its  anode. 
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